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1. Introduction 



Originally, the twisted mass formulation of lattice QCD |^ was invented to provide the Wil- 
son lattice Dirac operator with an infra-red cut-off for its eigenvalue spectrum in order to prevent 
the appearance of unphysically small eigenvalues when the quark masses become small. More re- 
cently, however, another quite astonishing property of twisted mass quarks has been discovered: in 
Ref. 1^ it was shown that pure Wilson twisted mass fermions lead to automatic G{a) improvement 
when tuned to full twist even without any improvement term added to the action. 

In simulations in the quenched approximation it has been demonstrated that the presence of the 
twisted mass term as an infra-red regulator is indeed a clear advantage when compared to the pure 
Wilson regularisation [Q]. However, this infra-red regulator appears, somewhat surprisingly, not to 
be crucial in simulations with light dynamical quarks. Here simulations even with standard Wilson 
fermions at pseudo-scalar masses below 300MeV are possible, at least when newer algorithmic 
techniques are used, see e.g. our own work of Ref. %. 

It appears therefore that for the case of dynamical simulations the aforementioned automatic 
6 [a) improvement of the (maximally twisted) Wilson formulation of lattice QCD becomes more 
important^ . In such a set-up only one parameter needs to be tuned in order to obtain &(a) improve- 
ment and there are no additional operator-specific improvement terms needed. Moreover, mixing 
patterns during renormalisation are greatly simplified. 

The disadvantage of the twisted mass formulation is the fact that the twisted mass term explic- 
itly breaks parity and flavour symmetry. The flavour symmetry breaking for instance introduces 
a mass splitting between the lightest charged and the uncharged pseudo-scalar mesons. However, 
all those effects stem only from the valence quark sector |^ (the sea quark determinant is flavour 
blind), are expected to be lattice artifacts of order &{a^). In particular, in a mixed action set-up, 
with for instance overlap valence fermions no breaking of parity or flavour chiral symmetry is 
expected. 

In the quenched approximation it is by now established that automatic ff{a) improvement 
works very well in practice ^, ^, [T^ . Also, the theoretical expectation of 0'{a^) lattice 



artifacts appearing in the pseudo-scalar mass splitting has been confirmed [|11|, [12p. In this pro- 
ceeding contribution we will provide first evidence that maximally twisted mass quarks are ^{a) 
improved also in dynamical simulations and that the residual 0'{a^) lattice artifacts are small. 
The lowest value of the pseudo-scalar mass covered in these simulations is significantly below 
300MeV, which is in a region where it is possible to confront the simulation data with chiral 
perturbation theory (;^PT) and determine therefore the validity range of ;tPT. Since simulations 
directly at the physical value of the pion mass are still somewhat too expensive, ;^PT could then be 
used to bridge the gap between the simulation data and the physical point. The low values of the 
pseudo-scalar mass reached in our work became possible only due to recent algorithmic develop- 
ments [|13|, ll, [T5|, |16[ I, l^, [18|, |19|] (See also the contribution by M. Clark to these proceedings). 



'Naturally, having an infra-red cut-off in the theory does guarantee a stable simulation, also in difficult situations. 
We consider therefore the presence of the twisted mass term still to be an advantage. In particular, when we think, e.g., 
of simulations much closer to the physical point than performed today, the presence of the twisted mass term may turn 
out to be crucial. 
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In particular, for all the Nf = 2 results reported in this proceeding contribution we have used our 
variant of the HMC algorithm described in Ref. |^. 

In addition to observables like pseudo-scalar mass and decay constant we will provide also 
first results for the aforementioned mass splitting in the pseudo-scalar sector and discuss the case 
ofNf = 2 + l + l flavours of quarks, i.e. taking the strange and the charm quark into account in the 
simulation. For recent summaries of the status of Wilson twisted mass fermions see Refs. [||,||], 



the review of Ref. [22] and references therein. 



2. Choice of Lattice Action 

The fermionic lattice action for two flavours of mass degenerate Wilson twisted mass quarks 
reads (in the so called twisted basis x) 

Stm = a'^Y.lxx [am + 4r + ijs -c^apClXx 

X I 

+ LZv Ux.v{r + yv)Xx+v + ul^yy{r-Yv)Xx-v ^ 
v=l -'J 



(2.1) 



where aniQ is the bare untwisted quark mass and a}X the bare twisted mass. T3 is the third Pauli 
matrix acting in flavour space and r is the Wilson parameter, which we set to one in our simulations. 

Twisted mass fermions are said to be at maximal twist if the bare untwisted mass is tuned to 
its critical value nicrit- We will discuss later on how this can be realised in practice. 

In the gauge sector we use the so called tree-level Symanzik improved gauge action (tlSym) 



[p3[], which includes besides the plaquette term Uu also rectangular (1x2) Wilson loops U^^\ 



Ss = ^Ubo i \{\-RtTv{Uu)}+bi i{l-ReTr([/i;2j} (2.2) 

with the bare inverse coupling j8, Z>i = — 1/12 and the normalisation condition Z^o = 1 — SZ^i. Note 
that with bi =0 this action corresponds to the usual Wilson plaquette action. 

Before presenting the actual results let us first briefly recall some earlier works which led to 
the present production set-up, see also Ref. [|^]. As mentioned before, in the quenched approxima- 
tion we have shown in an extended scaling test that ^(a) improvement works extremely well for 
maximally twisted mass quarks [Q, ^, ^]. In the context of this scaling test also several definitions 
for full twist were examined and it was found that in agreement with theoretical considerations 
[p4|, ESl E6|] the so called PC AC definition leads to small residual lattice artifacts of iffia^) in the 



whole range of - in particular, also small - masses investigated: for the PCAC definition, at a fixed 
value of ajJ. > 0, anient is determined from the condition that 

I,(3oAg(x)P-(0)) 

vanishes^ for large enough time separation such that a projection on the pion state is reached. This 
provides then anient as a function of ajj. and the final value of anient may be obtained by (linearly) 



^Note that mpcAC is sometimes denoted with mZ'^^'^, because the bilinears are defined in the twisted basis X- 
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L^xT 




K^rit(fl'/^min) 


ro/a(aAimin) 


3.9 


24^ X 48 


0.004 


0.160856 


5.184(41) 


4.05 


32^ X 64 


0.003 


0.15701 


6.525(101) 



Table 1: The parameters for the simulation. We denote by a/imin the smallest value of the twisted mass 
parameter ajJ. at which we have performed simulations. At this value of afi we determined the critical mass 
merit, or, equivalently the hopping parameter fQrit = 1 /(8 + lam^^iii). In order to convert to physical units, we 
take the value of ro/a as computed at a/Zmin- 



extrapolating amcrit(M) to = 0. However, the extrapolation of amcn\.{pi) to = is not necessary: 
taking the value of anicrit at the lowest value of a/i that is used in the simulation is an equally good 
definition. This remark is of particular importance for dynamical simulations where one might not 
be able to afford to determine amcnt at several values of ajx and then perform the extrapolation to 
At =0. 

Furthermore, in the dynamical case it is not clear at all how an extrapolation to a/x = could 
be performed: a generic property of Wilson type fermions is its non-trivial phase structure at finite 
lattice spacing with first order phase transitions close to the chiral limit, which was investigated 
in detail by means of twisted mass fermions in Refs. [27, 28, 29, 20] and which is in accordance 
with the picture developed in %¥! [3^ 31]. A very important consequence of this phenomenon 
is that at a non-vanishing lattice spacing the value of the pseudo-scalar mass cannot drop below a 
certain minimal value without being affected by the first order phase transition. From the analysis 
in xVT it is expected that the corresponding minimal value of the twisted mass parameter goes to 
zero proportional to towards the continuum limit. Hence, one can always find a largest value 
of the lattice spacing Omax where simulations with a given value of the pseudo-scalar mass can be 
performed. 

The actual value of ^max depends on many details, for instance on the gauge action used in 
the simulation. For example, when the Wilson plaquette gauge action is used one can estimate 
flmax ~ 0.07 — O.lfm to realise a pseudo-scalar mass of 250MeV. The situation improves, when 
also the 2x1 rectangular loops are included (see Eq. (2^)): the value of Omax increases when the 

The study of 



modulus of the coefficient b\ multiplying the rectangular part is increased [29 



different gauge actions led us to the conclusion that the tlSym gauge action, with a coefficient of 
^1 = —1/12, is a good choice to avoid on the one hand large effects related to the presence of the 
phase transition and on the other hand the possible conceptual issues concerning the DBW2 gauge 
action, with coefficient bi = -1.4088, pOl]. 



3. Set-up and First Results 

The main goal of the present work is to investigate, whether also in the case of Nf = 2 mass- 
degenerate flavours of dynamical quarks, the lattice artifacts are as small as has been found in the 
quenched approximation when maximally twisted Wilson fermions are used. The current status of 
our simulations are summarised in tables [l] and ^ While in table |T] we provide the parameters of 
our simulations, in table || we also give current estimates of the pseudo-scalar mass mps in physical 
units, where we used ro = 0.5 fm. In addition we provide the number of equilibration trajecto- 
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mps [MeV] 


^therm 


Mraj 




3.9 


0.0040 


280 


1500 


5000 


55(17) 




0.0064 


350 


1500 


5000 


23(07) 




0.0100 


430 


1500 


5000 


15(04) 




0.0150 


530 


1500 


5000 


06(02) 


4.05 


0.0030 


266 


1500 


2200 






0.0060 


372 


1500 


500 





Table 2: Current status of the simulations at )3 =3.9 and j3 = 4.05 using the simulation parameters listed in 
table |l|. In addition we give the current estimates of mps in physical units (ro = 0.5 fm) and of Tint(^') where 
available. 



ries A'^therm. the current statistics A'^traj we have available and estimates of the integrated plaquette 
autocorrelation time where available. 



The algorithm we used is a HMC algorithm with mass preconditioning [ ]13p and multiple time 
scale integration. It is described in detail in Ref. [j5|]. The trajectory length was set to 1/2 in all 
our runs. Our currently available estimates of the plaquette integrated autocorrelation time ZmtiP) 
quoted in table ^ are in units of T = 1/2. 

Full twist is realised in our simulations by tuning at a fixed value of j3 the PCAC quark mass 
"JpCAC to zero at the smallest value a/^min of the twisted mass parameter ajj,. Keeping then the 
value of Kent or, equivalently, merit, for all the other values of a/i provides ^(a) improvement for 
all relevant physical quantities (see Refs. [|25|, |8]]). We list in table [I] the values of Kent which 



we then kept fixed for all simulations at a given value of j8 . The dependence of ampcAC on ajj. is 
shown in figure [1(a). The fact that at values of ajj. > a/^min the PCAC quark mass ampcAC does not 



vanish wiU only lead to effects that are of ff{a ) and hence will not destroy i?{d) improvement. 

In order to make maximum use of the gauge configurations, we evaluate connected meson cor- 
relators using a stochastic method to include all spatial sources. This method involves a stochastic 
source (Z(2)-noise in both real and imaginary part) for all colours and spatial locations at one 
Euclidean time slice. By solving for the propagator from this source for each of the 4 spin compo- 
nents, we can construct zero-momentum meson correlators from any quark bilinear at the source 
and sink. This method involves 4 inversions per Euclidean time slice value selected since we chose 



to use only one stochastic sample. This "one-end" method is similar to that pioneered in Ref. [ |32| ] 



and implemented in Ref. [33]. We also employ a fuzzed source [|^] of extent 6 lattice spacings to 
enable study of non-local meson creation and destruction operators. 

In general, we save a gauge configuration every second trajectory and analyse meson corre- 
lators as described above from a selection of different Euclidean time slice sources. To reduce 
autocorrelations, we only use the same time slice source every 8-10 trajectories. Our primary sta- 
tistical error analysis is from a bootstrap analysis of blocks of configurations. This estimate was 
cross-checked by varying the block size. 

3.1 Setting the Scale 

As mentioned already we used for this first analysis the Sommer parameter ro with a value of 
0.5 fm to set the scale. While ro/a is measurable to high accuracy in lattice QCD simulations it 
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(a) ampcAC as a function of an al ji = 3.9. 



(b) ro as a function of afl al P = 3.9. 



Figure 1: PCAC quark mass and Sommer parameter as a function of the twisted mass. 

has the drawback that its value in physical units is not known very well, which might introduce a 
non-negligible uncertainty. Therefore, it could be advantageous to determine the scale using other 
quantities which are experimentally better known, such as fj^, /k or mx* which we plan to do in 
the future. 

In figure 1(b) we plot ro/a as a function of a jJ. at P = 3.9. Within the current errors the mass 



dependence of this quantity seems to be rather weak for the lowest three values ofajj. , which would 
allow for a constant extrapolation to the physical point (nijc = 137 MeV). Again, for the time being, 
we use the value of ro/a listed in table |l| determined at the value of a/imin as an estimate for ro/a 
at the physical point. 



3.2 /ps, mps and other observables 

In this subsection we present first results for /ps and mps- The pseudo-scalar mass nips is 
as usual extracted from the exponential decay of corresponding correlators constructed from local 
operators. In contrast to pure Wilson fermions, for maximally twisted mass fermions an exact 
lattice Ward identity allows to extract the pseudo-scalar decay constant /ps from the relation 



/ps = ^|(0|pi(0)|;r)h 



(3.1) 



with no need of computing any renormalization constant since it turns out that Zp = 1/Z^ [Q]. 
In order to account for finite size (FS) effects in mp$ and fp$ we used NLO ;tPT formulae as 



derived by Gasser and Leutwyler []35[]. The corrections decrease exponentially with a rate mps ■ L 
(which we take from our data) and are proportional to the overall factor (/ps -L)^^. The quantity 
/ps • L is estimated at the chiral point either via extrapolation of uncorrected data for /ps or from 
the phenomenological value of fj^ assuming ro = 0.5 fm. The difference between the results for the 
corrected /ps obtained in the two ways at the physical point is about 0.4 MeV, which is certainly 
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less than the uncertainty arising from neglecting higher orders in ;^PT for the FS effects. Using the 
socalled re-summed Luscher formulae with NNLO input (more precisely the data in Tables 



3 and 4 of Ref. [|36|]), we obtained for afys estimates of the finite size correction factors that are 
very similar to those coming from simple NLO ;|^PT, while for amps the Liischer formula leads to 
somewhat larger corrections. Note that for a definite statement about FS corrections, simulations 
with larger lattices would be needed in order to cross-check the predictions. Such simulations 
are in progress. 

In figure 2(a) we show (romps)^ (not FS corrected) as a function of the twisted mass ro/i for 
the two j8 -values available. The dependence is to a very good approximation linear in ro/x and the 
data at the two j3 values fall approximately on a common line. Note that the renormalization factor 
has not been taken into account here. The dotted line is not a fit to the data, but a line connecting 
the highest data point in rojj. of j8 = 3.9 with the origin and may serve to guide the eye. 

In figure 2{b] afps is plotted versus (amps)^ at j8 = 3.9 only. In the plot we give the FS 



corrected and uncorrected data. Clearly, the effect of the FS correction is visible but not huge. In 
any case, whether we perform a FS correction or not the data show a curvature. It is very tempting 
to perform a fit of the data according to the expectation from ^PT. For this, it would, however, be 
better to have a simulation at an additional value of ajj. > a/imin- Such a simulation is presently in 
progress and we will report the results in a forthcoming publication [p7|]. 

A linear extrapolation of the FS corrected data of /ps at j3 = 3.9 towards the physical point 
{m^ = 137 MeV) yields a value off,, = 129.0±0.5±1.0MeV (to be compared with/;^ = 131MeV). 
The first error comes from the errors on /ps, mps and the extrapolation in the quark mass, the sec- 
ond from the uncertainty on the value of r^/a. Systematic errors due to lattice artifacts, residual 
finite volume effects and the choice of ro = 0.5 fm are not included. 

In figure |3| we show a comparison of quenched and dynamical values of the average momentum 
(x) corresponding to the expectation value of a non-singlet operator in pion states. Although the 
data point from the dynamical simulation still has systematic uncertainties such as missing non- 
perturbative renormalisation and lack of a continuum extrapolation, the message of the graph is 
very encouraging: employing twisted mass fermions, it is clearly possible to address the approach 
to the physical point also for more complicated quantities than mps and fps - In addition it seems 
that such quantities can be computed rather precisely. 

3.3 Effects of Isospin Breaking 

As mentioned before, a very important issue to address is the effect of isospin breaking in the 
twisted mass formulation at finite lattice spacing. This effect is expected to be most significant 
in the mass splitting between the lightest charged and uncharged pseudo-scalar mesons. A first 
analysis at j8 = 3.9 and ajj. = 0.004, taking the disconnected contribution in the neutral channel 
fully into account, shows that the uncharged pseudo-scalar meson is about 25% lighter than the 
charged one. We get, in fact 

am^ =0.1361(5), a/7i° = 0.103(4) , 
or, expressed differently, r^{{m^)'^ — {m^Y) = c{a/ro)^ with c ^ —5, which has a magnitude that is 



a factor of 2 smaller than what was found in our quenched investigations [12]. Taking the value of 
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(a) (romps)^ as a function of ro/i at /3 = 3.9 (b) afps as a function of (amps)^ at j8 = 3.9 

and p = 4.05. without (squares) and with (circles) FS correc- 

tion. 

Figure 2: (romps)^ as a function of rojj. for both j3 values and a/ps as a function of (amps)^ for j3 = 3.9 
only. 



c = —5 and assuming a linear dependence of the charged and neutral pseudo-scalar mass-squared 
on ajJ., we can compute the value of ajJ. where the neutral pion becomes massless and the phase 
transition line starts. This critical value would be a/i^ ~ 0.0015 at j3 = 3.9 and hence well below 
the smallest value of ajj. employed in our simulations. If we consider the mass splitting in the 
vector meson sector, we find that this splitting is compatible with zero within our currently still 
large errors. Note that the uncharged pion being lighter than the charged one is compatible with 
predictions from ;^PT [^] in case the first order phase transition scenario is realised. 

Let us mention that the disconnected correlations needed for the meson are evaluated using 
a stochastic (Gaussian) volume source with 4 levels of hopping -parameter variance reduction [pT|]. 
We use 24 stochastic sources per gauge configuration and evaluate every 10 trajectories. 



3.4 Continuum Scaling of /ps 

In figure ^ we show the scaling behaviour of ro/ps with [a/rof using our current simulation 
data. These results shown in the plot are still affected by some uncertainties: First of all the values 
of rowjps are not (exactly) matched because there are only two values of apL available at j8 = 4.05. 
Second, we have taken - as mentioned before - the value of r^/a at a/^min at the two j3 -values as 
an estimate of the value at the physical point. 

However, we do not believe that those uncertainties are larger than the current statistical uncer- 
tainty. Therefore, figure Q suggests that lattice artifacts are rather small in ro/ps at the current level 
of accuracy. This result is a first indication that lattice artifacts are small and well under control in 
simulations with maximally twisted mass fermions, although further tests are clearly needed. 
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Figure 3: The average momentum of a parton in a pion. The blue squares correspond to quenched simula- 
tions using the Schrodinger functional. These simulations have controlled non-perturbative renormali- 
sation, finite volume effects and continuum limit. The black filled circles come from quenched calculations 
of the same quantity with maximally twisted Wilson fermions showing that with this approach much smaller 
quark masses can be reached, from Ref. [|9| (the open circle corresponds to a point where finite size effects 
may possibly still affect the data). The red filled circle represents a simulation point obtained with dynami- 
cal quarks at j3 = 3.9 (here the quenched renormalisation factor is used). Although a number of systematic 
effects, continuum limit and non-perturbative renormalisation still have to be analysed, it is to be noted that 
also at small quark masses dynamical results seem to agree with quenched ones. 



4. Including Strange and Charm Quarks 

Adding to a doublet of mass-degenerate twisted mass quarks the strange quark can be realised 
in two possible ways. Either, a single quark flavour is introduced, or, a non-degenerate doublet of 



quarks is taken, as suggested in Ref. [42]. We have followed the prescription of Ref. [42] where 
it has been shown that a real quark determinant with a mass-split doublet can be obtained if the 
mass splitting term is taken to be orthogonal in isospin space to the twist direction. The reason 
for choosing this solution is that we want to realise & (a) improvement also for the non-degenerate 
case by using maximally twisted Wilson fermions. In this way, again no further improvement 
coefficients have to be computed. 

The action for the non-degenerate quarks can be written in a form obtained from the physical 



basis of Ref. [42] by means of a chiral rotation around the direction 1 in isospin space to undo the 



twist of the Wilson term. The mass terms become 

Sna+i+i =a'^Xl{v^WN) + At5T^ + 'Mff75TVW} (4-1) 

X 

where mo has again to be set to its critical value. In this case the renormalised current mass of the 
two quarks is given by {pLc/Zp it /I5 /Z5), where Zp and Z5 are the renormalisation constans of the 
pseudo-scalar and scalar non-singlet operators (in the twisted basis). The i//^-fields couple to the 



gauge fields through the same Wilson-Dirac operator as the one of Eq. (2.1). 

In our Monte Carlo simulations we described the u and d quarks by a light mass-degenerate 
twisted mass doublet and the c and s quarks by a heavy non-degenerate twisted mass doublet. We 
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physical units at j3 = 3.9. 
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Figure 5: The squared pseudo-scalar pion and kaon masses as a function of the PCAC mass am^'^^'" of 
the light doublet (a/X/ = 0.0075, a^a = 0.2363, a/Xg = 0.2138, 16^ x 32 lattices with tlSym gauge action, 
j3 = 3.35). Owing to the presence of meta-stabilities, there is a gap between positive and negative values of 
am^'^^'^ . This implies a positive minimal value of the pion mass. 



used the PHMC updating algorithm [ |43| ] with multi-step polynomial approximation and stochastic 
correction in the update chain as described in Refs. [44, p^ ]. 

As has been discussed in Refs. ||6, 42, 46], in order to realise ^(a) improvement, tuning to full 
twist can be achieved by setting the values of jq. of the light and heavy doublet equal. Therefore 
with respect to the Nf = 2 simulations, we are left with only one additional parameter, i.e. the 
strange quark mass to be determined by the value of the kaon mass, while the charm quark is taken 



to be much heavier. In Ref. [ ]46| ] we give a detailed discussion of possible ways to determine the 
twist angle, renormalisation factors and meson masses. 

We find that simulations ofNf = 2 + l-\-l dynamical quarks are feasible in practice. The 
tuning effort is similar to the Nf = 2 case, modulo the need to know the ratio Zp/Zs to fix the 
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strange quark mass, and turns out to be possible without too much effort. Employing the tISym 
gauge action as before, we find for lattice spacings as coarse as a 0.15fm again signs of a first 
order phase transition which appears to be somewhat stronger than in the case of Nf = 2 degenerate 
quarks. As a consequence of this phase transition the minimal pion mass cannot be made to vanish. 
This is illustrated in fig. ^. Nevertheless, our simulations allow to estimate the values of j8 and /i 
where the lattice spacing is, say, a ~ O.lfm and safe simulations with a pion mass below 300 MeV 
become possible. ;^PT formulae for the case of non-degenerate quarks have already been developed 
[ p7| ] and have been confronted to our simulation results, see Ref. [ p6| ] for a discussion of this point. 

5. Discussion 

In this proceeding contribution we presented the status of an ongoing project to simulate two 
flavour lattice QCD with light maximally twisted mass quarks. Our first results are very encourag- 
ing. 

• We reach a pseudo-scalar mass of 280 MeV (setting the scale with ro = 0.5 fm) while the 
simulation is performing smoothly and stable. 

• Comparing results at 2 values of the lattice spacing, we find very good scaling, indicating 
that ^{a) improvement is at work and that also the ff{a^) lattice artifacts are small. This is 
in full accordance with earlier results in the quenched approximation. 

• We do see effects of isospin breaking which are most significant in the mass splitting of the 
neutral and charged pions and turn out to be 25%. This is already much smaller than the 
corresponding splitting obtained in the quenched approximation. Note that at our smaller 
value of the lattice spacing, corresponding to j3 = 4.05, the splitting would be at the 10% 
level only, if we assume an dependence of this effect as has been found in the quenched 
approximation. 

• We have performed first simulations with non-degenerate strange and charm quark flavours 
in addition to the mass degenerate light doublet. It turns out that tuning and, in particular 
reahsing maximal twist, in this case is not problematic and needs a similar effort as other 
approaches taking the strange quark into account. We have explored the phase structure 
in this case and determined those parameters where simulations would not be affected by 
unwanted (first order) phase transitions. Hence, calculations with twisted mass and Nf = 
2+1 + 1 flavours are well prepared. 

It has to be said that tuning to full twist has to be performed on the same large lattice as used later 
for the calculation of physical quantities. The reason for this is that we need to project to the pion 
state as fully as possible to determine the PCAC quark mass without being affected by finite size 
effects or excited state contributions. Thus the tuning itself is rather expensive. But, it only has 
to be done once, as it is the case for the determination of action improvement coefficients in other 
approaches. Note, however, that with twisted mass fermions we do not have to compute further 
operator specific improvement coefficients afterwards. 
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Our ongoing efforts concern the completion of the simulations at j3 =3.9 with an additional 
value of ajJ., which will be presented in a forthcoming publication [|7]|. The main thrust of the 
future work is to add another and finer value of the lattice spacing (and, maybe, also a coarser one). 
In this proceeding contribution we reported already about some first results in this direction which 
indicate only small scaling violations. Our large set of configurations is uploaded to the ILDG 
[ p8| , |49| ] and will be used for computations of more complicated physical quantities than discussed 
here. Another interesting direction we follow with the twisted mass approach is to use Osterwalder- 



Seiler or overlap fermions in the valence sector [50| and to perform simulations at non- vanishing 



temperature on which first results were presented at this conference [51 1. 
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